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Figure 1: Mesa, Shelf and 
Pedestal constraints are chosen 
to optimize an array factor with 
this sombrero shape. 


1 Introduction 

Suppose we have a planar transmit array of weighted antenna elements laid out on a lattice, 
a regular grid, and suppose each embedded element pattern is identical (typically achieved 
using guard elements). Then the far-held array pattern is just the product of the embedded 
element pattern and an array factor that is the 2D Fourier transform of the array taper. 
Mathematically, the array taper comprises impulses at the element locations with areas 
given by the element weights. We can speak of designing the weights, the taper, or the 
array factor. These are equivalent. 

Is it possible to optimize a transmit array taper that illuminates both a highly directive 
cluster of receive beams and a broad region of the sky, perhaps to support a ubiquitous radar 
function E El? Such a transmit array factor might allow us to track several targets and 
simultaneously monitor a broad region of interest for new targets. The goal of this report is 
to determine whether it is feasible to design a transmit array taper with such characteristics. 
The larger issues involved in designing such a radar system are outside the scope of this 
report and will not be addressed. 

The array factor that we desire has what we call a sombrero structure (Figure [T]), where 
the term mesa refers to the top of the sombrero, shelf refers to the top of the sombrero’s 
brim and pedestal refers to the bottom of the sombrero’s brim. The optimized amplitude- 
taper designs that follow have constraints in one or more of the angle regions defined by 
these terms. Constraints on the mesa region limit passband ripple, constraints on the shelf 
region limit the array-factor magnitude in the stopband, and constraints on the pedestal 
region define the height of the sombrero’s brim. 

In this report we are particularly interested in an “active” transmit array, in which each 
element is driven by its own power amplifier. There are two basic categories. In one cat¬ 
egory, the power radiated is limited by the power available from individual elements. The 
output-power capabilities of the individual amplifiers driving the elements is the determin¬ 
ing factor. In the other category, the power radiated is limited only by a constraint on total 
transmit power as determined by available prime power or by heat-dissipation issues, and 
individual amplifiers can provide power as needed. These two categories are idealizations; 
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in real systems power might well be constrained by both a total-power limit and by per- 
element power. However, consideration of these two extremes reveals a key choice: given 
a desired target illumination level, should our weight design minimize peak element power 
or total transmit power? 

Minimizing peak element power (without other design constraints) leads naturally to 
a phase-only taper. All element weights have the same magnitude. To minimize total 
transmit power, however, complex weights are needed, and design symmetries typically 
lead to weight phases being used only for beam steering so that the unsteered taper is a 
pure amplitude taper. 

Amplitude tapers are our main focus, not because they are what is most needed or 
desired but simply because we have first-rate tools for designing complex or amplitude 
tapers but are limited to the crudest of design methods for phase-only tapers. That situation 
is not entirely about software and methods but is about the problem itself: a phase is one 
degree of freedom, but a complex weight is two. Should we not expect that doubling the 
number of degrees of freedom will permit better optimization? Further and more important 
is that the two design spaces are different in a fundamental geometric sense. The amplitude- 
taper space is smoothly shaped and easily searched, while the phase-only space is jagged 
and rough, a problem only made incomprehensibly worse by the high dimensionality— 
often thousands—of the space. It’s largely unsearchable except by trial and error. 

So here we explore complex or amplitude tapers at some length and phase-only tapers 
quite tentatively. Phase-only tapers are a subset of complex tapers of course, so if there is 
no good solution to the complex-taper sombrero design problem, there cannot be a good 
solution to the phase-only sombrero problem either. Addressing complex tapers first then 
implicitly seeks to learn whether there is something in the nature of array factors themselves 
that precludes a sombrero structure. If our answer had been ”no sombreros allowed”, it 
would have applied to both types of tapers. As it turns out, that is not our answer. 

Certain details are common to every design presented in this report. All of the designs 
in this report assume a roughly circular 4507 element array arranged on a \/\f?> spaced 
equilateral triangular grid (Figure [2]). For ease of optimization, we require that the taper 
be unchanged when the whole array is subjected to any of the five indicated rotations or 
mirrored through any of the six lines shown in Figure[3} This automatically makes the array 
factor invariant to those same rotations and reflections. These array symmetries reduce the 
number of optimization variables by a factor of roughly 12. A real design might well 
call for a nonsymmetric pattern, but this approach provides the essential insight into the 
trade space at a much reduced computational cost. We assume a narrowband array and do 
not consider the wideband case. In each design we show boresight beams, however it is 
possible to phase steer any of these array factors off of boresight. 

Section [2] presents optimized amplitude taper designs. These amplitude taper designs 
start with a very simple pencil beam design and work towards the sombrero pattern de¬ 
scribed above. The report includes a progression of designs to show the effects of each 
constraint that is added to the optimization to eventually produce the sombrero pattern. 

Section [3]presents several phase-only tapers constructed as 2D chirps, with linear-chirp 
designs followed by nonlinear-chirp designs. Chirp parameters were chosen and tweaked 
by hand, simply by trial and error. The results are rough and serve mostly to illustrate the 
severity of the phase-only limitation. 
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Figure 2: Element layout 




Figure 3: The optimized tapers 
(and array factors) are made 
invariant to various reflections 
and rotations simply to make 
optimizing a large array man¬ 
ageable. 


2 Amplitude Tapers 

In a well-matched transmit array the power into each element is proportional to the squared 
magnitude of the corresponding element weight, so the total power input to the array is 
proportional to the sum of the squared weight magnitudes, the energy in the weight func¬ 
tion. In array optimization we are generally interested in maximizing the power density 
at some point in the far field, but should we do this assuming that the total power deliv¬ 
ered to the elements is fixed or that there is some limit to the power that can be applied to 
any individual element? The answer depends on the hardware design, and so each design 
was optimized both ways. For reasons of convenience we actually fix the far-field power 
density at some point in space and minimize one of these two measures of element power. 
Fixing the far-field power density and minimizing total/peak-element power is equivalent, 
to within a common scale factor on the weights, to fixing the total/peak-element power and 
maximizing the far-field power density. 

We use our MATLAB toolbox Opt OJ to design the amplitude tapers that follow. The 
Opt toolbox solves a particular class of convex-optimization problems in which a linear 
objective function is minimized subject to any number of constraints, each either a linear 
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constraint or a convex quadratic constraint. 

The first two designs in this section are reference cases, the Uniform design (page |4]) 
and the Mesa Reference design (page [8]). Following these reference cases are two series of 
optimized designs: designs minimizing the energy in the element-weight function ([Design 
|2.3| on page[l3lthrough |Design 2.9| on page[25]) and designs minimizing the peak magnitude 
of the element-weight function ( |Design 2.10| on page |27| through |Design 2.16] on page [39]). 
For each series of amplitude taper optimizations, we include 

• a simple, Taylor-like pencil beam design, 

• a design with mesa constraints, 

• a design combining the mesa constraints with a Taylor-like shelf constraint, 

• a design combining the Taylor-like pencil beam and the brim constraints, and 

• our final sombrero designs. 


Design 2.1 A simple reference case: test the code 

The reference design features uniform illumination. Each of the A weights is set to 1/A so 
that the array factor is unity, 0 dB, at boresight. No optimization is involved, but the code 
for computing and plotting the array factors and performance numbers is tested in a simple 
way here. 

This reference design serves a second purpose as well, one based on it possessing a very 
special property. This reference design can be represented by 2D discrete-time Fourier pair 
r n <-+ 12(f), where taper r n is a function of an integer two-vector element index n. In 
particular, r n = 1/A for n in the array and zero elsewhere. Using this reference design the 
signal power radiated in the boresight direction is proportional to |12(0) | 2 = 1, and the total 
power delivered to elements by the transmitter is proportional to ||r|| 2 = ]U n |r n | 2 = 1/A. 
Classically, boresight gain is taken as the ratio of these quantities |l?(0)| 2 /||r|| 2 = A. 
Likewise, if taper h n 21(f) is used instead, its gain relative to that reference gain is 


|ff(Q)l7llft|| 2 » 
I«(0)I7IMI 2 M 2 g 



n 


IWI 2 N 


I (ft. 71 2 

||/r|| 2 ||r || 2 


( 1 ) 


using on the right inner product (h, r) = h n r^. The point of this seemingly odd 
formulation of the gain improvement is that the ratio on the right cannot exceed unity, 
by the famous inequality of Cauchy, Schwarz, and Bunyakovsky. Therefore, no taper can 
improve on the boresight gain of this reference case, so for each design to follow we give 
the dB distance below this gain bound, classically referred to as taper loss, as 


(weight-energy taper loss) = — 101og 10 


\{h,r)\ 2 


( 2 ) 





















This reference case also serves as the gain bound in the peak-element power sense. If 
we replace the element energy ||/i|| 2 with the peak element magnitude // max | 2 , 


|tf(o)| 2 /|/w | 2 

|i?(0)| 2 /|r max | 2 


n 

\K^N~ 2 ' 


(3) 


But the numerator is equal to (h, r) 2 /|r max | 2 , so for each design we can report the dB 
distance below the gain bound as 


(max-weight taper loss) 


—101og 10 


l(MI 2 

|/imax| 2 iV 2 |r m ax | 2 


(4) 


Array factors are plotted in what is often and somewhat confusingly called either sine 
space or direction-cosine space. The far-held plane wave radiated with 3D radian vector 
wavenumber 27rk has amplitude proportional to <f>(k)f7((kA — s) T B), where <f>(k) is the 
embedded element pattern and H ((kA — s) T B) is the phase steered array factor of interest. 
Here 77(f) is the 2D discrete-time Fourier transform of the element weights h n , s is a unit 
vector in the steering direction, and B is a dimensionless 3x2 matrix characterizing the 
geometry of the element-layout grid. For a plane wave, normalized wavenumber kA is a 
unit vector, but difference kA — s is not and can have a vector length approaching two. For 
this reason, we plot 77 (k T AB) without restriction to the unit-vector case. Decomposition 
k = k + k^ into array plane and perpendicular boresight components allows us to take 
advantage of the orthogonality of k ± to the columns of B, both of which always lie in 
the array plane, through the resulting fact that k T B = k T B. To plot 201og 10 |77(k T B)| 
then, we will plot 20 log 10 |77(k T B) | versus components of k in the azimuth and elevation 
directions, and most often we will write the array factor as |77(k T B) | as well. 

The array factor is displayed at three different zoom levels in the first column of each 
figure. 


galaxy zoom: A view of the array factor showing several periods. 


world zoom: A big-picture view of the array factor. The latitude and longitude lines repre¬ 
sent elevation and azimuth respectively, and the interior of the “world” is the visible 
region, the part of the array factor corresponding to physical directions. Elevation 
and azimuth lines are spaced at 15° intervals. 


nation zoom 60 dB: A close-up view of the array factor showing details of the transition 
region. The plot spans —7.5° to 7.5° in elevation and azimuth, with a color scale of 
60 dB. Elevation and azimuth lines are spaced at 1° intervals. The dB color scale uses 
contrasting colors for the — 3 dB and — 2 dB levels in order to make the associated 
beamwidths visible. 


nation zoom 4 dB: A close-up view of the array factor showing details of the mesa, with 
a color scale of 5 dB. The dB color scale uses contrasting colors for the — 3 dB and 
—2 dB levels in order to make the associated beamwidths visible. 
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slice plot: An array-factor slice plotted on both a linear and dB scale. The slice shown is 
from boresight to the edge of the visible region along elevation. This plot depicts 
qualitative behavior of the array factor and is particularly helpful when the optimiza¬ 
tion behavior is unexpected. 

element weight vs. radius: The element weights are plotted versus the radius. Due to 
the symmetry of the problem specifications, none of the optimal amplitude weights 
have nonzero imaginary components, so all phases are 0° or 180°. Negative and 
nonnegative element weights are different colors (or will be for designs with weights 
of both polarities). 

element weights: The element weights are plotted in 3D perspective. (This plot is not 
particularly revealing for this design, but will be for others to follow.) 

The reported taper loss and taper efficiency numbers are two different ways of expressing 
the same measurement. They are not two different measurements. 
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Results: 


weight-energy taper loss (weight-energy taper efficiency) 0 dB (100%) 
max-weight taper loss (max-weight taper efficiency) 0 dB (100%) 


Figure 4: Uniform element weights 
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Design 2.2 A mesa reference case 

The taper losses of (|2j) and (|4]) jump upward when we add a mesa to the designs because of 
the large increase in noise-gain factor 


E 


\h n \ 2 = j \H(f)\ 2 df = A 2 |B t B| 1/2 f \H(k T XB)\ 2 dk, 

J unit J one 


'unit 

square 


one 
period 


where di is differential area in two-vector f and dk is differential area in the array-plane 
component k of three-vector k. The integral on the right becomes much larger when 
H( k T AB) is given a near-unity value over the entire mesa rather than over a pencil beam 
as before. What this means is that the |Design 2.1 [ reference is a poor choice for tapers with 
mesas: How do you know if a taper loss of 12 dB is good or not? It turns out that we can 
extend the notion of taper loss to explicitly account for a mesa in a way that is particularly 
clean for amplitude tapers. 

How do we characterize the mesa mathematically? Array factors are plotted as func¬ 
tions of the array-plane component k of three-vector k, so the mesa is actually a set of 
array-plane values. Each such value can be written as a weighted combination of a pair of 
linearly independent array-plane basis vectors. The columns of basis matrix B are such a 
pair, but it is convenient here to instead use the rows of B + , the Moore-Penrose pseudoin¬ 
verse of B that is formally given as B + = (B T B) -1 B T and that is realized in matlab with the 
pinv () function. Its rows are in the array plane because each is just a linear combination 
of the basis-vector columns of B, and this new basis is appealing because if we represent 
an array-plane value k using a row two-vector f of coordinates so that k T = yfB + , the 
simplifying relationship k T AB = fB + B = f results. 

Choosing the mesa now amounts to choosing a set of f vectors. We will soon be tak¬ 
ing averages across the mesa, and averages are nicely handled by the expected values of 
probability theory, so let us suppose that f, and implicitly k as well, is a random vector dis¬ 
tributed probabilistically across one array-factor period and, in particular, with probability 
concentrated in the mesa. Then 


average 
amplitude gain 


£|tf(k T AB)| E\H(i)\ 


By distributing probability in the mesa according to the weighting desired in the average, 
we could in principle obtain a weighted average. We do not do so, however. 

Consider some specific cases. 


Case: single-point mesa. If f = 0 with probability one, 


average 
amplitude gain 


|ff(0)l 

IWI 


For this degenerate mesa the average amplitude gain is just the boresight gain used to 
derive weight-energy taper loss in |Design 2. 1| (which begins on page [4]). There is no 
reason to compute average amplitude gain separately with the single point considered 
as a mesa. 
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Case: mesa defined by a pdf. This is the main case of interest. Suppose the probability 
of f falling into a subset T of the unit square is 

P[feT] =J P (f)df, 

where df is differential area in the f plane so that p(f) on the unit square is a proba¬ 
bility density function (pdf) for f. Then let p({) be periodic with the unit square as 
its period, so that it can be inverse Fourier transformed. After inverse transforming 
\H (f) | as well, we can use 2D Fourier pairs h n <-+ \H(f ) | and </>„ p( f) along with 
Parseval’s relation and \\h\\ 2 = \h\ | 2 , which is obvious in the frequency domain, to 
write 

amplitude gain = pjj jj H ( f )I P( f ) di ( 5 ) 

square 

= (h,_<f>) 

\M ' 


There are two subcases. 

Subcase: zero-phase taper. The amplitude tapers of this report not only give each 
weight in the taper a zero phase, the design symmetries imposed also give each 
point of the associated array factor a zero phase. Though somewhat confusing, 
the custom in signal processing is to refer to tapers with the latter property as 
zero-phase tapers (or filters) or, less precisely, as linear-phase tapers (or filters). 
Typically a zero-phase taper not only has //(f) real everywhere in the mesa, 
generally in the mesa it also has //(f) > 0 so that //(f) | = //(f ) and 

average _ 1 f „/ f x / f x , f _ (h, (fi) 

amplitude gain “ 11 h 11 JJ? ,rjp,rj " r ~ ||/,|| ' 

square 

Since h n = 0 for n not corresponding to actual array elements, inner product 
(h , 0) = has many zero terms. We can just as well decompose 0 as 

(fin = &n + (fi±n, where 0 n and (fi_ n can be nonzero only for n corresponding to 
elements actually inside and outside the array respectively, and write ( h , <fi) = 
(h, (fi). It is then useful to consider the ratio 

average 

amplitude gain _ (h, cfi) 

m = wM’ 

because by the Cauchy-Schwarz-Bunyakovsky inequality the quantity on the 
right will never exceed one, but we can make it equal to one by setting h = cfi. 
We characterize the closeness of a mesa design to this optimum, using termi¬ 
nology loosely (because technically we are averaging the gain and not the loss), 
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by reporting 


average 


(average mesa taper loss) 



( 6 ) 



with confidence that in this zero-phase subcase it must be nonnegative. 

Subcase: a taper with phase variation. Later, in Section[3](beginning on page|4T|) on 
phase-only tapers, we will continue to report (| 6 j) but with the average amplitude 
gain calculated by numerically evaluating the integral in (|5]). In this subcase we 
have no Cauchy-Schwarz-Bunyakovsky bound on the values this average mesa 
taper loss might take. 

The Airy taper. We give our zero-phase tapers mesas by defining pdf p(i) = p(k T AB) 
to be constant inside a circle in k of mesa radius Ay and to be zero outside that circle. An 
inverse 2D discrete-time Fourier transform of p( f) and subsequent restriction of (f) n to the 
array aperture then gives (j) n as a sampled version 


<£ n = 2 jinc(27rAy||BAn||) 


(7) 


of the Airy pattern from optics, using a jinc function defined using a Bessel function: 

for x 7 ^ 0 , 



for x = 0 . 


The scaling has been chosen so that y 0 = / square p(f) di = 1 as required. 

Given that our array size is fixed, Ay is the one free parameter in the definition of o n . 
Average mesa taper loss is presented for every taper in this report for which a mesa is an 
intended design feature, and for all amplitude tapers a common mesa radius parameter Ay is 
used in calculating <£ n . Just as we do not expect a filter’s ripple bandwidth, 3 dB bandwidth, 
and noise bandwidth to be the same, here we do not assume that Ay should exactly match 
the design radius used to specify the mesa size in the optimization constraints. Instead we 
determine this common value of Ay by computing the average mesa taper loss of [Design 
2 . 4 | for a range of Ay values. The result is shown in Figure[5j The lowest average mesa taper 
loss is obtained with Ay = 0.0882/A, so that value is used for calculating average mesa 
taper loss for each amplitude-taper design that has a mesa. 

Of course taper loss is indifferent to scaling, so rather than take the scaling implied 
by h = cj), we can if we wish set the boresight gain to approximately unity. The scaling 
of Fourier pair (f) n ■<—> p( f) yields f sqmre p( f) di = 1, so if we scale (J) n , d n , and p(k T AB) 
by |B T B| 1 // 2 7 r(AyA) 2 , after scaling the last of the three will have unit height in the mesa. 
After scaling, the first of the three becomes our reference taper and will transform to an 
array factor of approximately unit mesa height but featuring truncation-related ripples and 
sidelobes. 
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Figure 5: Average mesa ta¬ 
per loss of |Design 2Aj versus 
ip r , the angular radius of the 
mesa. We choose the direction- 
cosine-space mesa radius k r = 
sin(^v)/A in 0 to minimize 
this |Design 2~4| average mesa 
taper loss. 
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Results: 


weight-energy taper loss (weight-energy taper efficiency) 13.80 dB 
max-weight taper loss (max-weight taper efficiency) 29.09 dB 
average mesa taper loss (average mesa taper efficiency) 0 dB 


(4.17%) 

( 0 . 12 %) 

( 100 %) 


Figure 6: Reference design for all amplitude tapers with mesas 
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Design 2.3 Weight-energy minimization, Taylor-like pencil beam 


The first set of optimizations minimize the array-factor energy (or weight energy, thanks to 
Parseval). Note that for the examples with mesas this is not exactly the same as minimizing 
the weight-energy taper loss since the boresight gain is not explicitly fixed. Nor is it equiva¬ 
lent to minimizing average mesa taper loss, because the average mesa gain is not explicitly 
fixed. This can lead to some slightly unexpected results when comparing designs, as shown 
later in the discussions for |Design 2.8| and |Design 2.12[ However, the peak constraints on 
the mesa do approximately fix both the boresight and average mesa gains, and thus the net 
result is to approximately minimize the two taper losses that depend on total weight energy. 

This design is a pencil beam with a Taylor-like response. The following constraints are 
applied to the filter (• on the array-factor plots): 


• upperbound the array-factor sidelobes by —35 dB over the shelf region, from 2.5° 
outward. 

• constrain the array factor at boresight to 0 dB. 

The element weights are nonnegative. 

The points where constraints are applied are visible on each of the array-factor plots (•). 
Due to the symmetries in the problem formulation (Figure[3]), these points, which are spread 
across roughly 1/12 of the array-factor period, constrain nearly the entire universe. Of 
course if we were to remove those required array-factor symmetries from the formulation, 
for example to enable the inclusion of some asymmetric sidelobe limits, the number of 
variables optimized would be much larger. 
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Results: 


weight-energy taper loss (weight-energy taper efficiency) 
max-weight taper loss (max-weight taper efficiency) 

weight-energy level 


0.83 dB 
6.25 dB 
-35.71 dB 


(82.62%) 

(23.72%) 


Figure 7: Weight-energy minimization of Taylor-like pencil beam 
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Design 2.4 Weight-energy minimization, mesa constraint 

Next, we test the mesa constraint on our design. The mainbeam is broadened to 4.25° so 
that a cluster of approximately 19 pencil beams will fit inside the mesa (Figure [8]). The 
following constraint is applied to the filter (•) on the array-factor plots): 

• bound the array factor between — 1 dB and 0 dB in the mesa region, from boresight 
to 4.25°. 


The weight-energy taper loss is 14.99 dB. This large number results from the dissim¬ 
ilarity between the present wide beam and the reference pencil beam of |Design 2.1| The 
average mesa taper loss is a more appropriate measure because it compares to the mesa- 
reference wide beam of |Design 2.2| The average mesa taper loss is 0.02 dB. That it is 
very small is no surprise, because in Figure [5] we fixed the mesa radius of |Design 2.2| to 
minimize the weight-energy taper loss of this very design. The element weights are both 
negative and nonnegative. 

If desired, one could eliminate the hexagonal rays that extend from the mainbeam by 
modifying the shape of the circular array. 



Figure 8:19 beams packed into 
a hexagonal cluster. 
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Results: 


weight-energy taper loss (weight-energy taper efficiency) 
max-weight taper loss (max-weight taper efficiency) 
average mesa taper loss (average mesa taper efficiency) 

weight-energy level 


14.99 dB 
29.97 dB 
0.02 dB 
-22.55 dB 


(3.17%) 

( 0 . 10 %) 

(99.47%) 


Figure 9: Weight-energy minimization with a mesa constraint 
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Design 2.5 Weight-energy minimization, mesa constraint, —35 dB shelf 
constraint 

Here we combine the mesa constraint with the Taylor-like shelf constraint to get a broader 
beam with a Taylor-like response. The following constraints are applied to the filter (•) on 
the array-factor plots): 

• bound the array factor between — 1 dB and 0 dB in the mesa region, from boresight 
to 4.25°. 

• upperbound the array-factor sidelobes by —35 dB over the shelf region, from 6.8° 
outward. 
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Results: 


weight-energy taper loss (weight-energy taper efficiency) 
max-weight taper loss (max-weight taper efficiency) 
average mesa taper loss (average mesa taper efficiency) 

weight-energy level 


15.11 dB 
30.98 dB 
0.06 dB 
-22.43 dB 


(3.08%) 

(0.08%) 

(98.60%) 


Figure 10: Weight-energy minimization with a mesa constraint and a —35 dB Taylor-like 
shelf constraint 
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Design 2.6 Weight-energy minimization, mesa constraint, —23 dB shelf 
constraint 

Here, we change the shelf height from — 35 dB, a reasonable shelf height for a Taylor-like 
array, to —23 dB, the target height of the brim for our sombrero designs. The following 
constraints are applied to the filter (•) on the array-factor plots): 

• bound the array factor between — 1 dB and 0 dB in the mesa region, from boresight 
to 4.25°. 

• upperbound the array-factor sidelobes by —23 dB over the shelf region, from 6.8° 
outward. 
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Results: 


weight-energy taper loss (weight-energy taper efficiency) 
max-weight taper loss (max-weight taper efficiency) 
average mesa taper loss (average mesa taper efficiency) 

weight-energy level 


15.00 dB 
30.20 dB 
0.02 dB 
-22.54 dB 


(3.16%) 

( 0 . 10 %) 

(99.52%) 


Figure 11: Weight-energy minimization with a mesa constraint and a —23 dB Taylor-like 
shelf constraint 
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Design 2.7 Weight-energy minimization, Taylor-like pencil beam, 30° 
brim constraints 

Here we go back to a pencil beam design (we remove the mesa constraints) and add the 
brim constraints (the pedestal and shelf constraints in Figure [TJ to produce a sombrero¬ 
shaped array factor with a pencil-beam mesa. The following constraints are applied to the 
filter (•) on the array-factor plots): 

• constrain the array factor at boresight to 0 dB. 

• upperbound the array-factor sidelobes by —23 dB over the shelf region, from 2.5° 
outward (•) on the array-factor plots). 

• lowerbound the array factor by —25 dB in the pedestal region, from 6.8° to 30° (o on 
the array-factor plots). 
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Results: 


weight-energy taper loss (weight-energy taper efficiency) 
max-weight taper loss (max-weight taper efficiency) 

weight-energy level 


8.13 dB 
37.01 dB 
-28.41 dB 


(15.37%) 

( 0 . 02 %) 


Figure 12: Weight-energy minimization of a Taylor-like pencil beam with 30° brim con¬ 
straints. 
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Design 2.8 Weight-energy minimization, sombrero constraints, 30 ° brim 


Next, mesa constraints are reintroduced to generate our desired sombrero pattern. The 
constraints for the sombrero pattern are: 

• bound the array factor between — 1 dB and 0 dB in the mesa region, from boresight 
to 4.25° (•) on the array-factor plots). 

• upperbound the array-factor sidelobes by —23 dB over the shelf region, from 6.8° 
outward (•) on the array-factor plots). 


• lowerbound the array factor by —25 dB in the pedestal region, from 6.8° to 30° (o on 
the array-factor plots). 


The constraints of this sombrero design include all of the constraints from the Taylor- 
like Mesa design ( |Design 2T6| ) plus one additional (pedestal) constraint. Therefore the set 
of possible (feasible) optimization outcomes of this sombrero design is a subset of the set 
of possible optimization outcomes of the Taylor-like mesa design. One would expect 

quality (best((small set))) < quality (best((large set))) 

if (small set) C (large set). It appears superficially that whether this expectation is met 
depends on the quality measure chosen. The 

, . (l5.11dB Taylor-like Mesa design, 

(weight-energy taper 1oss)<ib — j 1493dB Sonlbrcro design . 


and 


(energy ) dB = 


—22.43 dB Taylor-like Mesa design, 
—21.66 dB Sombrero design. 


The weight-energy taper loss results do not match our expectations, while the weight- 
energy levels do meet our expectations. The energy is actually minimized by the opti¬ 
mization, so the results there are not surprising. Relationship ([!]) can be written in decibel 
form as weight-energy taper loss to energy is 


(weight-energy taper loss) dB = 

(ideal gain) dB + (energy) dB — 20 log 10 1 (boresight array factor) |. (8) 

Ideal gain is not affected by the optimization, but the boresight-array-factor term can be. 
Here 

,,. . . _ ,, ( 1.00 dB Taylor-like Mesa design, 

-20 logjo | (boresight array factor) 1= | 0 .05 dB Sombrero design. 

If we use this to adjust the taper loss. 


(weight-energy taper loss) dB 

— 20 log 10 1 (boresight array factor) | 


14.00 dB Taylor-like Mesa design, 
14.88 dB Sombrero design, 


which behaves as expected. 

We conclude from this array-factor design that it is possible to optimize a sombrero 
pattern for an array. 
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Results: 


weight-energy taper loss (weight-energy taper efficiency) 
max-weight taper loss (max-weight taper efficiency) 
average mesa taper loss (average mesa taper efficiency) 

weight-energy level 


15.86 dB 
40.55 dB 
0.75 dB 
-21.68 dB 


(2.60%) 

( 0 . 01 %) 

(84.23%) 


Figure 13: Weight-energy minimization with sombrero constraints (30° brim) 
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Design 2.9 Weight-energy minimization, sombrero constraints, 45 ° brim 

Here is another example of the sombrero pattern, with the brim extending out to 45°. The 
constraints for the sombrero pattern are: 


• bound the array factor between — 1 dB and 0 dB in the mesa region, from boresight 
to 4.25° (•) on the array-factor plots). 

• upperbound the array-factor sidelobes by —23 dB over the shelf region, from 6.8° 
outward (•) on the array-factor plots). 

• lowerbound the array factor by —25 dB in the pedestal region, from 6.8° to 45° (o on 
the array-factor plots). 
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Results: 


weight-energy taper loss (weight-energy taper efficiency) 
max-weight taper loss (max-weight taper efficiency) 
average mesa taper loss (average mesa taper efficiency) 

weight-energy level 


16.32 dB 
44.77 dB 
1.20 dB 
-21.22 dB 


(2.33%) 

(0.003%) 

(75.86%) 


Figure 14: Weight-energy minimization with sombrero constraints (45° brim) 
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Design 2.10 Max-weight optimization, Taylor-like pencil beam 

The second set of optimizations minimize the peak magnitude of the element weights. 
From ([3]) we can obtain a relationship analogous to ([8]) but for the max-weight taper loss: 


(max-weight taper loss) dB = 

(ideal gain) dB + 10 log 10 (N \h max | 2 ) - 20 log 10 1 (boresight array factor) |. (9) 

This design is a pencil beam with a Taylor-like response. The following constraints are 
applied to the filter (•) on the array-factor plots): 

• upperbound the array-factor sidelobes by —35 dB over the shelf region, from 2.5° 
outward. 


• constrain the array factor at boresight to 0 dB. 

Comparing to the weight-energy optimization with identical constraints ( |Design 2~3j ), 
note that the peak optimization buys us nothing (there is no visible change in the peak 
coefficient) and costs us weight energy. 
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Results: 


weight-energy taper loss (weight-energy taper efficiency) 2.16 dB (60.75%) 

max-weight taper loss (max-weight taper efficiency) 5.32 dB (29.35%) 

max-weight level —67.75 dB 


Figure 15: Max-weight optimization of a Taylor-like pencil beam 
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Design 2.11 Max-weight optimization, mesa constraint 

Next, we test the mesa constraint on our design. The mesa constraints broaden the main- 
beam to 4.25° to fit approximately 19 pencil beams inside the mesa (Figure[8]). The follow¬ 
ing constraint is applied to the filter (•) on the array-factor plots): 

• bound the array factor between — 1 dB and 0 dB in the mesa region, from boresight 
to 4.25°. 

The average mesa loss of this case is significantly worse than the complimentary weight- 
energy optimization case ( |Design 2~4| ). Although the max-weight taper loss is reduced, we 
also have very large and undesirable sidelobes surrounding the mesa. 
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Results: 


weight-energy taper loss (weight-energy taper efficiency) 
max-weight taper loss (max-weight taper efficiency) 
average mesa taper loss (average mesa taper efficiency) 

max-weight level 


22.58 dB 
22.62 dB 
8.07 dB 
-51.46 dB 


(0.55%) 

(0.55%) 

(15.58%) 


Figure 16: Max-weight optimization with a mesa constraint 
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Design 2.12 Max-weight optimization, mesa constraint, —35 dB shelf 
constraint 

Here we combine the mesa constraint with the Taylor-like shelf constraint to get a broader 
beam with a Taylor-like response. The following constraints are applied to the filter (•) on 
the array-factor plots): 

• bound the array factor between — 1 dB and 0 dB in the mesa region, from boresight 
to 4.25°. 

• upperbound the array-factor sidelobes by —35 dB over the shelf region, from 6.8° 
outward. 


Comparing this design to pesign 2,5[ we find the surprising result that both the weight- 
energy and max-weight taper losses have improved, albeit marginally. This simply reflects 
that here the array factor peaks up at boresight, while in |Design 23] the array factor is 
depressed at boresight. 
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Results: 


weight-energy taper loss (weight-energy taper efficiency) 
max-weight taper loss (max-weight taper efficiency) 
average mesa taper loss (average mesa taper efficiency) 

max-weight level 


14.37 dB 
28.99 dB 
0.13 dB 
-44.14 dB 


(3.66%) 

(0.13%) 

(97.05%) 


Figure 17: Max-weight optimization with a mesa constraint and a —35 dB Taylor-like shelf 
constraint 
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Design 2.13 Max-weight optimization, mesa constraint, — 23 dB shelf 

constraint 

Here, we change the shelf height from — 35 dB, a reasonable shelf height for a Taylor-like 
array, to —23 dB, the target height of the brim for our sombrero designs. The following 
constraints are applied to the filter (•) on the array-factor plots): 

• bound the array factor between — 1 dB and 0 dB in the mesa region, from boresight 
to 4.25°. 

• upperbound the array-factor sidelobes by —23 dB over the shelf region, from 6.8° 
outward. 


Here again the improvement in max-weight taper loss over pesign 2.6| is minimal and 
accompanied by a significant loss in the average mesa taper loss. 
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Results: 


weight-energy taper loss (weight-energy taper efficiency) 
max-weight taper loss (max-weight taper efficiency) 
average mesa taper loss (average mesa taper efficiency) 

max-weight level 


16.20 dB 
28.01 dB 
1.26 dB 
-46.07 dB 


(2.40%) 

(0.16%) 

(74.76%) 


Figure 18: Max-weight optimization with a mesa constraint and a — 23 dB Taylor-like shelf 
constraint 
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Design 2.14 Max-weight optimization, Taylor-like pencil beam, 30 ° brim 
constraint 


Here we go back to a pencil-beam design (we remove the mesa constraints) and add the 
brim constraints (the pedestal and shelf constraints in Figure [TJ to produce a sombrero¬ 
shaped array factor with a pencil-beam mesa. The following constraints are applied to the 
filter (•) on the array-factor plots): 

• constrain the array factor at boresight to 0 dB. 

• upperbound the array-factor sidelobes by —23 dB over the shelf region, from 2.5° 
outward (•) on the array-factor plots). 

• lowerbound the array factor by —25 dB in the pedestal region, from 6.8° to 30° (o on 
the array-factor plots). 
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Results: 


weight-energy taper loss (weight-energy taper efficiency) 11.45 dB 
max-weight taper loss (max-weight taper efficiency) 30.22 dB 

max-weight level —42.86 dB 


(7.17%) 

( 0 . 10 %) 


Figure 19: Max-weight optimization of Taylor-like pencil beam with 30° brim constraints 
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Design 2.15 Max-weight optimization, sombrero constraints, 30 ° brim 

Next, mesa constraints are reintroduced to generate our desired sombrero pattern. The 
constraints for the sombrero pattern are: 


• bound the array factor between — 1 dB and 0 dB in the mesa region, from boresight 
to 4.25° (•) on the array-factor plots). 

• upperbound the array-factor sidelobes by —23 dB over the shelf region, from 6.8° 
outward (•) on the array-factor plots). 

• lowerbound the array factor by —25 dB in the pedestal region, from 6.8° to 30° (o on 
the array-factor plots). 

The improvement in max-weight taper loss over pesign 2.8| is modest. 
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Results: 


weight-energy taper loss (weight-energy taper efficiency) 
max-weight taper loss (max-weight taper efficiency) 
average mesa taper loss (average mesa taper efficiency) 

max-weight level 


16.21 dB 
35.12 dB 
1.96 dB 
-37.96 dB 


(2.39%) 

(0.03%) 

(63.70%) 


Figure 20: Max-weight optimization with sombrero constraints (30° brim) 
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Design 2.16 Max-weight optimization, sombrero constraints, 45 ° brim 

Here is another example of the sombrero pattern, with the brim extending out to 45°. The 
constraints for the sombrero pattern are: 


• bound the array factor between — 1 dB and 0 dB in the mesa region, from boresight 
to 4.25° (•) on the array-factor plots). 

• upperbound the array-factor sidelobes by —23 dB over the shelf region, from 6.8° 
outward (•) on the array-factor plots). 

• lowerbound the array factor by —25 dB in the pedestal region, from 6.8° to 45° (o on 
the array-factor plots). 


As in many of the previous designs, compared to pesign 2.9| we see a small improve¬ 
ment in the max-weight taper loss and similar increases in the weight-energy and average- 
mesa taper losses. 
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Results: 


weight-energy taper loss (weight-energy taper efficiency) 
max-weight taper loss (max-weight taper efficiency) 
average mesa taper loss (average mesa taper efficiency) 

max-weight level 


16.46 dB 
38.65 dB 
2.20 dB 
-34.43 dB 


(2.26%) 

( 0 . 01 %) 

(60.29%) 


Figure 21: Max-weight optimization with sombrero constraints (45° brim) 
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3 Phase-Only Tapers 

This section is based on 0| but includes both more examples and material on taper losses 
that is not present there. 

Earlier taper designs in this report all featured real element weights. Optimization of 
real weights or even complex weights is practical even for large problems because the 
design constraints are mathematically convex: given any two tapers u n and v n that satisfy 
the constraints, any taper obtained by linearly interpolating between them, like .7 u n + .3 v n 
or .lw n + .9 v n or .4 u n + .6 v n , satisfies those constraints as well. For such a convex 
optimization problem, any local optimum is also the global optimum. Further, there are 
efficient algorithms for finding that global optimum. 

However, robust and reliable optimization of a phase-only array taper is practical only in 
the most limited of special cases, simply because a constant-amplitude constraint \h n \ = 1 
itself is nonconvex. Interpolation between two unit-magnitude complex numbers does not 
in general yield a unit-magnitude result; a straight line between two points on the circum¬ 
ference of a circle does not remain on that circumference. A large array will have several 
thousand of these nonconvex constraints, one per element except when elements are related 
by symmetry. For large arrays this typically gives the resulting optimization problem a vast, 
even ghastly number of local optima, and there is no reliable way of discovering which of 
them is the global optimum. 

Viable phase-only strategies then tend to either use a heuristic method to pick a rea¬ 
sonable suboptimal design or to use a preliminary design chosen in that way to initialize 
a follow-on optimal search for a nearby local optimum. Examples of heuristic design ap¬ 
proaches can be found in (51© 7, [8| ( 9 , [TOj. Local optimization approaches are presented in 
mmmmmmm. Nonlocal optimization methods, such as simulated annealing 
and genetic algorithms, have also been used although they do not guarantee global 
solutions either. There are also hybrid approaches such as llT9ll . in which a small number 
of phase basis functions is heuristically chosen and the corresponding coefficients are then 
locally optimized. Whether any of these procedures yields a useful result depends very 
much on the particulars. 

A simple FM-chirp heuristic approach. One heuristic approach commonly used for 
uniform line arrays in ID borrows from classical FM waveform design. The core idea is 
simply to hope that the magnitude of a chirp’s Fourier spectrum at a given frequency / 
is proportional to the amount of time the instantaneous frequency, given for chirp e j<9 ^ 
by /inst(f) = ^^ 7 ^, spends near /. It’s very much like an integral change of variable 

/ = /inst(f) or t — /i~t(/) so that dt = df. The chirp spends time dt between 

instantaneous frequencies / and / + df, so we might hope that its Fourier spectrum will 

l rlf — 1 ( f \ j 

have magnitude at least roughly proportional to | |. The messy part mathematically 

is generally dealing with the inverse-function aspect of this, but there are classic cases in 
which the relationships are not difficult to work out, and exploring them brings out the 
limitations of the approach. 
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For example, for t € (— the FM chirp waveform 


u>(£) 


e M/»)(t/r) f()r | t | < t / 2 , 

0 otherwise. 


has instantaneous frequency f3t/T, which sweeps at a constant rate across interval (—f, f) 
as t increases. One might hope then that its Fourier magnitude spectrum is roughly constant 
on that interval and zero elsewhere. In fact it is approximately so when time-bandwidth 
product /3T is large and only then. In that case, 


mn\ 


some constant for |/| <(3/2, 
0 otherwise. 


Use unsampled 2D tapers and their transforms for simplicity. We begin the extension 
to 2D by establishing a simplified notation. Relating that new notation to our existing 
notation is a tad tedious but is along the lines of the usual development of sampling in 
signals-and-system texts. Let the orthonormal columns of 3 x 2 matrix A point in the 
azimuth and elevation directions when looking out along the array boresight vector. If two- 
vector x contains azimuth and elevation coordinates, we can express a position in the array 
plane as Aa?, and we can express the array-plane component k of wavenumber k using a 
two-vector k of azimuth-elevation coordinates as k = A k. 

For simplicity let’s then assume the taper is obtained by sampling a continuous func¬ 
tion (j{x) at the values of x corresponding to element positions. The sampling process 
can be represented in the x and k domains and continuous-”time” (really space) Fourier 
transformation in 2D, where each sum is over all two-vectors of integers, as 

9(x) 2D ^ T G(k) 

x * 

A 2 |B T B|^<5(a; - A T ABn) 2 ^ T ^S(k - ±(A T B)- T m) 

n m 

ll 

A 2 |B T B|^(A T ABn) 5(x - A T ABn) 2D ^ T e~ i2nkTAJxBn 

n V ^ n 

h n ^ y 

tf(fe T A T AB) =tf(k T AB). 

The transform in the last line follows by direct evaluation of the transform integral. The 
underbraces identify the actual taper (left) and array factor (right) in terms of the 2D 
discrete-time (space) Fourier pair h n H( f) used earlier for amplitude tapers. Here 
g(x ) must have units of 1/(length) 2 to make h n dimensionless, and the actual array factor 
f/(k T AB) = H( k T AB) is then related to the Gik) that we will use here by k = A k and 

tf(fc T A T AB) = ^ G(k — i(A T B)- T m). 

m 

We actually have no further need of such details of scaling and periodic replication, as 
we derive results only up to a scale factor and assume always that the sampling-related 
replication amounts to little more than periodic extension. 
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Brick-wall beams from phase bowls. Extending the ID chirp idea to a rectangle in 
array-plane coordinates is now straightforward. Beamwidth and time-extent matrices (3 
and T are diagonal with positive elements on those diagonals: 


g( x ) 


1 e j7TX T (3T 1 x 
some area 

o 


for T 1 x G S, 
otherwise. 


( 10 ) 


Here set S is a 1 x 1 square centered on the origin. In 2D the instantaneous (really 
“instanspaceous” since x comprises position coordinates, but we lack such a word) fre¬ 
quency of a signal e) 0lx> is just k m st = ^VO(x), with the gradient taken with respect 
to coordinate vector x. The phase of signal g(x) of (jTOj) has the parabolic “bowl” shape 
O(x) = Trx J f3T~ 1 x and instantaneous frequency fc msl = /3T -1 a?, and the latter ranges 
across those fc msl for which /3 -1 fei nst G <S. And indeed, when matrix (3 T has large diagonal 
elements, 


\G(k)\ 


some constant for (3 1 k G S, 
0 otherwise. 


This continuous phase-only taper has a rectangular outline and gives an approximately 
rectangular beam in spatial frequency k (sine space). 

For an outline elliptical in vector x we could use 


g( x ) 


1 £j7ra? T /3T x x 
some area 

0 


for x T T 2 x < |, 
otherwise, 


where 2x2 symmetric time and beamwidth matrices T and f3 here are positive definite. 
The major and minor diameters of the taper are given by the eigenvalues of T. Again 
fcjnst = (3T~ l x, and in the ellipse where g(x) is nonzero it ranges over the ellipse in fc inst 
for which kj nst f3~ 2 ki nst < 1/4, so 


|G(*)I 


some constant 
0 


for k 1 (3 2 k < 
otherwise. 


resulting in an elliptical beam with major and minor diameters given by the eigenvalues of 
(3. Other support regions can be used, resulting in beams with matching geometries. 

Any of these continuous-space functions can be sampled spatially to obtain element 
weights for an array. The primary problem with such “FM beams” is that the approxima¬ 
tion error is a function of the (eigenvalues of the) space-beamwidth product /3T, which 
(eigenvalues) must (all) be large. Even the largest radar arrays do not have enough ele¬ 
ments for this approximation to be reasonable—it takes hundreds of thousand elements to 
get down to a few dB of ripple—so we are reduced to tweaking parameters when (each 
eigenvalue of) /3T is rather small in hopes of “lucking into” a decent array factor. 


Design examples approximating brick-wall beams. We illustrate in the design exam¬ 
ples to follow by presenting several actual array factors. Our array and desired beam here 
are both circular, so the elliptical case above applies. That circularity makes both (3 and 
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T -1 into scaled identity matrices, so we can parameterize their product as oil = /3T _1 
using a single scalar parameter a. Each of our first group of phase-only tapers therefore 
takes the form 

e j7ra||a; 


Our spatial-extent matrix T is fixed by the array geometry, so choosing our one available 
parameter a corresponds to choosing beamwidth matrix (3. Given a desired beam size, 
we have no degrees of design freedom left. Here we simply pick several representative a 
values for the sake of illustration. 

We do notice that the boresight value of the array factor seems especially sensitive to a, 
so given a general range for a that gives us roughly the beamwidth we seek, we generally 
end up fine tuning a to set this boresight value to something nonextreme. 


for 11 a? 11 in our 4507-element array, 
otherwise, 


( 11 ) 



A crude attempt at a Gaussian beam. The phase-bowl design examples just discussed 
have instantaneous frequencies fcj nst that are linear in position x and so implicitly aim at 
having their array factors approximate “brick wall” spectra. It is no surprise therefore that 
nasty array-factor ripples arise as a result, and it certainly seems reasonable to attempt some 
sort of gradual rolloff as an alternative, in hopes that the approximations will be better, that 
ripple magnitudes will be significantly lower. 

Towards this end, consider the taper 


9(x) 



for ||a;|| < r 0 , 

otherwise. 


( 12 ) 


Different authors define the error function with different scalings, but here z = erf 1 (y) 
refers to the odd function defined for positive arguments by 

y = erf(z) = 4y [ e~ u2 du. 

Jo 

A little algebra shows instantaneous frequency as a function of x to be 

fcinst(*) = 2 ^Vphase(*) = |j^* y^erf _1 (||*||/r 0 ), 

and change of variable k = k mst (x) then results (shown in 0]) in the ratio of differential 


areas 


Then 


= 2r o ( ^ erf(||fe||/y^fcg) \ e _|| fc ||2 /(2fco2) 
dk 7 rk%\ WkWUPiM ) 


7 


(13) 


|G(fc)| » (some constant) x ^ erf(W ^> e -ll*HV(«g), (14 ) 

at least for sufficiently high r 0 k 0 , which here serves the role of a time-bandwidth product. 
The quantity in large parenthesis in ( |T3| ) goes to unity as k goes to zero, and the error 
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function approaches for large arguments, so the Fourier spectrum approximated in 

( |T4| ) is essentially Gaussian for \\k\\ <C k 0 but rolls off somewhat faster than Gaussian for 

||fc||»fc 0 . 

We present designs below for r 0 = 20.31A, and several different values of k 0 . It is no 
surprise that the array factor approximates the value here more closely when r 0 k 0 is large 
and that for beam sizes of interest r 0 k 0 in fact is fairly small. Still the resulting array- 
factor ripple magnitudes can sometimes be made far lower than for the phase-bowl tapers 
by choosing r 0 and k 0 carefully. Generally approximation quality is far more sensitive to k 0 
than to r 0 , and again the value of the array factor at boresight, at the origin, is particularly 
sensitive. Very small deviations from the k 0 values shown can triple or quadruple mainbeam 
ripple. 


Sombrero-shaped array factor. One possible way to generate a sombrero-shaped phase- 
only array taper is to replace ||*|| in ( [12] ) with some function /(||*||) giving 


g( x ) 



for ||a;|| < r 0 , 

otherwise. 


(15) 


We present several designs below using /(||*||) = b\\x\\ + (o m + ||*|| m ) 1 / m — a. As with 
the previous phase-only designs, small deviations in any one of the five parameters k 0 , r 0 , 
a, b, or m cause major changes in the array factor. 


Presentation of individual designs. For each phase-only design that follows, the array- 
factor magnitudes are shown normalized to a maximum magnitude of 0 dB. Array-factor 
phases are nonzero but are not shown. The phase of the element weights are plotted in 
3D in the lower right. The plot above that represents instantaneous frequency along an 
arbitrary cut through the array origin with k- UKl equal to the plotted quantity times a unit 
vector in the cut direction. The average mesa taper loss of ((6]) (page [TO]) is calculated using 
<£ of {7]) (page [TO]) with a mesa radius parameter k r customized to that specific design by 
searching for a minimum in the spirit of Figure |5] (page ]TT|). 
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Design 3.1 Linear-FM chirp, small a 

Here we use LFM chirp ( [IT] ) with a small a. The 3 dB beamwidth of this design is similar to 
the beamwidth of our amplitude-taper pencil-beam designs ( |Design 23] and [Design 2.10| ), 
but the rolloff in this design is much slower. 
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Parameters: 


a = 0.003/A 2 
k r = 0.0705/A 





Results: 

weight-energy taper loss (weight-energy taper efficiency) 6.45 dB (22.64%) 

max-weight taper loss (max-weight taper efficiency) 6.45 dB (22.64%) 

average mesa taper loss (average mesa taper efficiency) 1.11 dB (77.36%) 

Figure 22: Linear-FM chirp, small a 
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Design 3.2 Linear-FM chirp, medium a 

Here we use LFM chirp ( fTT| ) with a medium a. The 3 dB beamwidth of this design is similar 
to the beamwidth of our amplitude-taper mesa constrained designs ( [Design 24| and [Design 
|2.11| ). Compared to those mesa constrained designs, the ripple in this design is predictably 
bad, the rolloff is much slower and the overall sidelobe level is higher. However, we do see 
a significant improvement in the max-weight taper loss when compared to our optimized 
amplitude-taper designs. 
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Results: 


weight-energy taper loss (weight-energy taper efficiency) 
max-weight taper loss (max-weight taper efficiency) 
average mesa taper loss (average mesa taper efficiency) 


15.04 dB (3.13%) 
15.04 dB (3.13%) 
0.84 dB (82.45%) 


Figure 23: Linear-FM chirp, medium a 
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Design 3.3 Linear-FM chirp, large a 

Here we use LFM chirp < [TT| ) with a large a, which gives us a design with a larger mesa. 
Both the ripple in the mesa and the difference in the boresight array factor are very sensitive 
to a. 
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Results: 

weight-energy taper loss (weight-energy taper efficiency) 29.46 dB (0.11%) 
max-weight taper loss (max-weight taper efficiency) 29.46 dB (0.11%) 
average mesa taper loss (average mesa taper efficiency) 0.58 dB (87.59%) 

Figure 24: Linear-FM chirp, large a 


51 





















































































































































Design 3.4 Nonlinear-FM chirp for narrow Gaussian beam 

Here we use ( fl2l ) with a fixed r 0 and a small k 0 . 
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Results: 

weight-energy taper loss (weight-energy taper efficiency) 11.15 dB (7.68%) 

max-weight taper loss (max-weight taper efficiency) 11.15 dB (7.68%) 

average mesa taper loss (average mesa taper efficiency) 1.25 dB (74.98%) 

Figure 25: Nonlinear-FM chirp for Gaussian beam, small fc 0 
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Design 3.5 Nonlinear-FM chirp for medium-width Gaussian beam 

Here we use ( fl2l ) with a fixed r 0 and a medium k 0 . This design is comparable to [Design 
|3.2| and our amplitude-taper mesa designs ( [Design 2~4] and |Design 24~Tj ). Compared to 
Design 2.4| the weight-energy and average-mesa taper losses are both somewhat worse, and 
as expected the max-weight taper loss is significantly improved. Compared to pesign 2.11 
all the taper losses are significantly improved. The ripple in the mesa is much smaller in 
this design than in|Design 3.2|and the roll-off and sidelobe levels are similar. 
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r 0 = 20.31A 
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:• jj: 


weight-energy taper loss (weight-energy taper efficiency) 16.04 dB (2.49%) 
max-weight taper loss (max-weight taper efficiency) 16.04 dB (2.49%) 
average mesa taper loss (average mesa taper efficiency) 1.07 dB (78.14%) 


Figure 26: Nonlinear-FM chirp for Gaussian beam, medium k 0 
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Design 3.6 Nonlinear-FM chirp for wide Gaussian beam 

Here we use ( fl2l ) with a fixed r 0 and a large k 0 . These parameters generate a design with 
a larger mesa comparable to |Design 3J| The ripple in the mesa is much smaller in this 
design than in pesign 3.3[ but the roll-off and sidelobe levels are not as good. 
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Results: 

weight-energy taper loss (weight-energy taper efficiency) 24.94 dB (0.32%) 
max-weight taper loss (max-weight taper efficiency) 24.94 dB (0.32%) 
average mesa taper loss (average mesa taper efficiency) 0.84 dB (82.46%) 

Figure 27: Nonlinear-FM chirp for Gaussian beam, large k Q 
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Design 3.7 Nonlinear-FM chirp for sombrero beam, first try 

Here we use <tnj). For this design, the mesa is flat for approximately 2.5°, there is a —15 dB 
shelf to approximately 30° and a sidelobe level of —30 dB at 90°. 
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Results: 

weight-energy taper loss (weight-energy taper efficiency) 18.26 dB (1.49%) 
max-weight taper loss (max-weight taper efficiency) 18.26 dB (1.49%) 
average mesa taper loss (average mesa taper efficiency) 1.32 dB (73.85%) 

Figure 28: Nonlinear-FM chirp for sombrero beam, first try 
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Design 3.8 Nonlinear-FM chirp for sombrero beam, better transition 

Here is another design using ( fl5] ). 
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Parameters: 


k 0 = 0.3291/A 
r 0 = 20.31A 
a = 10A 
6 = 0.4 
m = 6 

k r = 0.4352/A 




Results: 


weight-energy taper loss (weight-energy taper efficiency) 
max-weight taper loss (max-weight taper efficiency) 
average mesa taper loss (average mesa taper efficiency) 


15.32 dB (2.94%) 
15.32 dB (2.94%) 
1.37 dB (72.95%) 


Figure 29: Nonlinear-FM chirp for sombrero beam, better transition 
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Design 3.9 Nonlinear-FM chirp for sombrero beam, pointy beam 

Here we use This design has a pointy beamed mesa with a slow transition to a shelf 
and rolloff comparable to|Design 3.7 
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Results: 


weight-energy taper loss (weight-energy taper efficiency) 
max-weight taper loss (max-weight taper efficiency) 
average mesa taper loss (average mesa taper efficiency) 


15.78 dB (2.64%) 
15.78 dB (2.64%) 
1.35 dB (73.36%) 


Figure 30: Nonlinear-FM chirp for sombrero beam, pointy beam 
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4 Conclusions 


Optimizing sombrero-shaped amplitude (or complex) tapers is, although not a trivial prob¬ 
lem, one that we can successfully solve with our current optimization tools. Of the two 
amplitude-taper optimization objectives presented, minimizing total transmitter power while 
(approximately) fixing boresight or mesa-averaged radiated power density generally pro¬ 
duces better results than minimizing the maximum element power. For the most part the 
latter simply increases the weight-energy and average mesa taper losses without bringing 
the max-weight taper loss anywhere near a reasonable number and at the additional cost 
of irregular and sometimes large sidelobes. There may be specific situations where some 
gentle constraints on the peak element weights would be useful, but in general if the array 
to be designed is heavily peak limited then an amplitude taper appears to be a rather poor 
choice. Unfortunately, current digital-array driver concepts are inherently peak limited. 

Phase-only tapers do, when designs comparable to optimized amplitude tapers are dis¬ 
covered, significantly improve the taper loss in the peak-limited element case. However, we 
cannot generate a sombrero shaped phase-only design with the same degree of control and 
performance as an amplitude taper when the same array geometry is used. The phase-only 
sombrero shaped designs shown here were generated by tweaking five different parameters 
in ( p~5] ). Each of these five parameters has significant effects on the array factor, but those 
effects are not systematic. The designs were a result of guessing a starting point and hand 
adjusting the parameters to attempt to design a reasonable sombrero shaped array factor. 
This process could be formalized by directly optimizing such parameters or even the phases 
themselves, but the presence of large numbers of local minima would make it unpractically 
expensive to exhaustively search for the global minimum. 
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Design 

Optimization 

variable level 

weight- 

energy 

taper 

loss 

max- 

weight 

taper 

loss 

average 

mesa 

taper 

loss 

Uniform 

N/A 

N/A 

0 dB 

0 dB 

N/A 

Taylor Pencil 

weight energy 

-35.71 dB 

0.83 dB 

6.25 dB 

N/A 

Mesa 

weight energy 

-22.55 dB 

14.99 dB 

29.97 dB 

0.02 dB 

Taylor Mesa —35 dB 

weight energy 

-22.43 dB 

15.11 dB 

30.98 dB 

0.06 dB 

Taylor Mesa —23 dB 

weight energy 

-22.54 dB 

15.00 dB 

30.20 dB 

0.02 dB 

Taylor Brim 

weight energy 

-28.41 dB 

8.13 dB 

37.01 dB 

N/A 

Sombrero 30° 

weight energy 

-21.68 dB 

15.86 dB 

40.55 dB 

0.75 dB 

Sombrero 45° 

weight energy 

-21.22 dB 

16.32 dB 

44.77 dB 

1.20 dB 

Taylor Pencil 

max weight 

-67.75 dB 

2.16 dB 

5.32 dB 

N/A 

Mesa 

max weight 

-51.46 dB 

22.58 dB 

22.62 dB 

8.07 dB 

Taylor Mesa —35 dB 

max weight 

-44.14 dB 

14.37 dB 

28.99 dB 

0.13 dB 

Taylor Mesa —23 dB 

max weight 

-46.07 dB 

16.20 dB 

28.01 dB 

1.26 dB 

Taylor Brim 

max weight 

-42.86 dB 

11.45 dB 

30.22 dB 

N/A 

Sombrero 30° 

max weight 

-37.96 dB 

16.21 dB 

35.12 dB 

1.96 dB 

Sombrero 45° 

max weight 

-34.43 dB 

16.46 dB 

38.65 dB 

2.20 dB 


Table 1: Amplitude-taper results 
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weight- 

max- 

average 



energy 

weight 

mesa 



taper 

taper 

taper 

Design 

variable 

loss 

loss 

loss 

Linear-FM chirp 

a = 0.003/A 2 

6.45 dB 

6.45 dB 

1.11 dB 

Linear-FM chirp 

a = 0.006/A 2 

15.04 dB 

15.04 dB 

0.84 dB 

Linear-FM chirp 

a = 0.015/A 2 

29.46 dB 

29.46 dB 

0.58 dB 

Nonlinear-FM chirp 
Gaussian beam 

ko = 0.07/A 

11.15 dB 

11.15 dB 

1.25 dB 

Nonlinear-FM chirp 
Gaussian beam 

k 0 = 0.1068/A 

16.04 dB 

16.04 dB 

1.07 dB 

Nonlinear-FM chirp 
Gaussian beam 

k 0 = 0.3291/A 

24.94 dB 

24.94 dB 

0.84 dB 


k 0 = 0.3291/A 




Nonlinear-FM chirp 
Sombrero beam 

ro = 24A 
a = 12.048A 
b = 0.697 
m = 5 

18.26 dB 

18.26 dB 

1.32 dB 


ko = 0.3291/A 




Nonlinear-FM chirp 
Sombrero beam 

r 0 = 20.31A 
a = 10A 

6 = 0.4 
m = 6 

15.32 dB 

15.32 dB 

1.37 dB 


ko = 0.3291/A 




Nonlinear-FM chirp 
Sombrero beam 

ro = 24A 
a = 12A 
b = 0.68 

15.78 dB 

15.78 dB 

1.35 dB 


m = 6 


Table 2: Phase-only taper results 
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